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Abstract 

An integro-differential equation on a tree graph is used to model the 
evolution and spatial distribution of a population of organisms in a river 
network. Individual organisms become mobile at a constant rate, and 
disperse according to an advection-diffusion process with coefficients that 
are constant on the edges of the graph. Appropriate boundary condi- 
tions are imposed at the outlet and upstream nodes of the river network. 
The local rates of population growth/decay and that by which the organ- 
isms become mobile, are assumed constant in time and space. Imminent 
extinction of the population is understood as the situation whereby the 
zero solution to the integro-differential equation is stable. Lower and up- 
per bounds for the eigenvalues of the dispersion operator, and related 
Sturm-Liouville problems are found, and therefore sufficient conditions 
for imminent extinction are given in terms of the physical variables of the 
problem. 

1 Introduction 

The problem of persistence of a population of organisms in an environment with 
predom i nantly unidirectional flow gives rise to the so-called "drift paradox" 
( Miilleil . [l98i . For riverine habitats, the problem is to determine physical, bio- 



logical or dynamical mechanisms by which a species whose organisms spend time 
drifting in the water column, avoids being driven to extinction. Recently, a body 
of qu antitative work has been used to address this problem, s ee ([Lutscher et al . 



2005h and references therein. In particular in lLutscher et all (|2005f ) the authors 



used the following integro-differential equation to model spatio-temporal dy- 
namics of the population, at low-density values, on a river stretch of length 



du 



— {x,t)^ru{x,t)-^iu{x,t) + ^ij IC{y,x)u{y,t)dy. (1) 
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Here, u{x,t) is the number of individuals per unit length at point x and time 
r > is a net population growth rate at small densities, /i > is the rate 
at which individuals become mobile, and K. is a. "dispersion kernel", namely 
/C(y, x) is the probability that a mobile individual disperses from y to x. 

Of main interest is to find conditions under which the trivial solution u = is 
a stable state of equation ^ . Under those conditions we say that the population 
faces "imminent extinction" , namely it cannot endure low population density 
values. The complementary situation, where s tability of t he zer o solution to ([1]) 



does not hold, is referred to as "persistence" (|j0rgenseij . 120041 ). 

The integro-differential Equation ([1]) is of the Barbashin type, and a neces- 
sary condition for the stability of the trivi al solution is that its largest Lyapunov 



exponent be negative I Appell et at |2000| ) . Namely, 



r - + ^sup{Re(w); w G Spec(/C)} < 0, (2) 

where the spectrum Spec(/C) of the integral operator /C[/] — )C{y, x)f{y) dy 
denotes the collection of all its eigenvalues, that is, complex numbers oj such 
that JC[f] — ujf for some / in the appropr i ate fu nction space. 

The analytical work in Lutscher et al ( 20051 ) yielded the following results: 



(i) If the dispersal kernel does not depend on the size of the environment I, 
the largest eigenvalue ujic of /C is an increasing function of I, and therefore, 
there exists a "critical domain size" I under which inequality ([2|) does not hold, 
and the population faces imminent extinction, (ii) If the motion of mobile 
individuals is assumed to follow an advection-diffusion process with drift v and 
diffusion coefficient D > 0, then upstream dispersal is likely enough to ensure 
the existence of persistence scenarios, even for high values of v. These results 
rely on the assumption that the environment is a single river stretch and ignore 
boundaries of the environment; in particular the conditions in (ii) are found for 
an advection-diffusion kernel on the whole real line trunca t ed to [0, /]. 

In this paper we extend the analysis of iLutscher et al (l2005l) in two direc- 



tions. First the finiteness of the environment is taken into account and particular 
boundary conditions are considered. Secondly, equation ([T]) is considered in a 
directed tree graph made of several connected segments with different phys- 
ical properties, and hence arrive at conditions for persistence of populations 
undergoing advection-diffusion in a river network. The goal is to lay the foun- 
dation for a deeper understanding of the role that some of the properties of real 
river networks, e.g. heterogeneity, connectedness, scaling, etcetera, play on the 
existence of persistence conditions. 

The rest of the introduction is devoted to revisiting the conditions of stability 
on the one-dimensional case, and summarizing the results in the case of a net- 
work. In section [5] the equations and operators for advection-diffusion on river 
networks are derived. Section ([3]) deals with the relationship between JC and 
a related Sturm-Liouville operator, along with some important consequences. 
Lastly, sections dU and ([5]) are devoted to proving the stability conditions for 
populations on networks. 
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1.1 The one-dimensional case 



As a motivation, and to fix ideas, assume the population density evolves on a 
river stretch [0,1] as in ([T]). Consider also the following particular dispersion 
mechanism. Let P(y,x,t) be the transition probabilities of a diffusion process 
with constant diffusion D, drift v; with absorbing boundary condition at a; = 0, 
and reflecting boundary conditions at x ~ I. Then P satisfies the backwards 
equation, 

dP d^P dP , , dP , , 

Suppose now that mobile organisms follow the diffusion process described in 
for an exponentially distributed random time of mean — units of the time scale 
t defined by ([T]) . Then, the dispersion kernel is 



IC{y,x) 



/ ae''^*P{y,x,t)dt. (4) 
Jo 



Estimates on the eigenvalues of K, can be given in terms of the following two 
non-dimensional quantities: 

^ vl ^ V 

The first one being the classical Peclet number, measuring the relative impor- 
tance of advection and diffusion in the dispersion process. The variable Q com- 
pares the mean velocity of the channel with the velocity needed to transverse 
the whole channel during the mean time individuals are mobile. 

The results in the one-dimensional case follow from standard theory. How- 
ever, details are included here to illustrate some of the ideas that will be later 
used in the case or river networks. 

Theorem 1.1. The largest eigenvalue is given by ujjc = where vi satisfies 

l + iQP<.i<l + iQP + ^^. (6) 

Proof. Any eigenvalue for the operator /C in ([4]) is given by w = -i where v is & 
g-eigenvalue of the associated Sturm-Liouville operator 

an - i-pf')' + qf, P{^) ■■= q{x) ^pix). (7) 

Namely, ly satisfies C[u] — vqu for some u G C'^([0, 1]) such that m(0) = u'{l) — 0. 
In particular, u is has the form u{x] v) — _Be"^ -1- Ce^^ for _B, C G M, and a, (i 
given by 



V + Jv"^ - ADaiv - 1) ^ V - Jv"^ - ADa{v - 1) 
2D ' ^'■= 2b ■ 
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Evaluation of the boundary conditions yields that both a and /? must have 
non-zero imaginary part, and so the estimate on the left hand side of ([6]) holds. 
Moreover, i^i can be obtained as the smallest solution to 

tan(/6H) + ^ 0, ^ = ^ ^ ^p^Q ' + 1 (8) 



where b{i') = 2I) ~ 4Z?ct(j^ — 1)|. The first of the equations on ([5]) defines 
6 as a function of I, and one can differentiate implicitly to obtain -gj < 0. Also 
d| ^ 2^ ^ follows that < 0. The upper bound in © follows simply 

from noting that since ^ > 0, the smallest positive solution to tan(^&) + ^ = 
must satisiy lb G {^,n). □ 

The following corollaries contain the extinction and critical domain condi- 
tions that follow from theorem 

Corollary 1.2. 1. If r > fi the population will persist. 

2. Given v, D, a and I, there is a critical population growth rate rcrit '■= 
1 — such that if r < rcrit, the population will face imminent extinction. 
Moreover, 

PQ rent 4Q 
4+PQ /i P(Q2 +7r2) +4Q' ^ ^ 

3. Given < r < /i and fixing all other parameters, there is a critical length 
Icrit such that if I < Icrit, the population will face imminent extinction. 

Figures ([1]) and ^ show the dependence of the critical values in terms of 

2 

non-dimensional quantities: P, Q, PQ — jj^ and ^. 

1.2 Dispersion on networks. 

By a network T we understand a directed, finite, binary, geometric graph em- 
bedded in IR2 assume that each edge e £ T allows a sufficiently smooth 
parametrization, contains no self-intersections, and is finite, therefore can be 
considered as the interval e — [0,le\- A point on F is then denoted as the pair 
(e, x) with ^ X ^ le. The "root edge" is denoted as r, and each edge has either 
zero or two children edges connected at the point (e,Ze). At each endpoint of 
an edge is located a node of F. The set of nodes is A^(F) and boldface is used 
to denote individual nodes. The "upstream node" of edge e is e = {e,le)', the 
"root node" of F is </> = (r,0). 

The children edges of r are (0) and (1). Inductively, the children edges of 
some edge e = (i) are (iO) and (il), also denoted as eO and el. The upstream 
node of e has three possible representations 

e = (e,/e) = (eO,0) ^ (el,0). 
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log(P) 

Figure 1: Values of rcrit as a function of P and Q. Note that the horizontal 
scale is logarithmic. The dashed line shows the boundary of positive critical 
reproductive rate. 
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Figure 2: Critical Peclet number Pdit- Given - and fixed v, D and a, the value 

Pcrit = "^jj"* gives the length of channel below which extinction is imminent. 
Note the horizontal logarithmic scale. The dashed line corresponds to the lower 
bound on the left hand side of (0) 



6 



We say that e is a "leaf edge" if it has no children edges. Upstream nodes of 
leaf edges form the upstream boundary of F, denoted as U{T). The "boundary" 
of r is defined as dV = U{T) U {4>}. Internal nodes are J(r) = N{T) \ dV. 

Values of a function / : F — M are denoted by fe{x) = /(e, x). That is, fe is 
the restriction of / to the edge e. Conversely, given any collection of functions 
{fe : [0,/e] ^ M; e e r} with the property Uih) = /eo(0) = /ei(0) for all 
e G /(r), one can construct the function / = J^eer /e^e on T. The derivative 
of a function / : F K at a point {e,x) in the interior of e is understood 
in the usual sense. For edge endpoints (e, 0) and (e, le) the derivative /' is 
understood as the right and left derivatives respectively. At an internal node e, 
the derivative of a function is not well defined as it can take three values. For 
nodes e S dT, /'(e) denotes the appropriate one-sided derivative. 

The space of functions that are continuous inside each edge of F and at all 
nodes, is denoted by C(f ). The set C(F) contains, in contrast, those functions 
that are continuous for (0, ^e) for each edge, but such that their values at nodes 
may not correspond to the one-sided limit taken along an incoming edge. Sim- 
ilarly, C"(F), n = 1, 2, . . . , oo denotes the space of functions that are n times 
continuously differentiable in the interior of all edges. 

A tree graph F will serve as a model for a river network with outlet located 
at <p, and each edge is the one-dimensional representation of a stream. Given 
an edge e, the following strictly positive variables are defined; water velocity Ve, 
cross-sectional area A^, and diffusion coefficient D^- We are then assuming that 
these correspond to representative spatio-temporal mean values of the physical 
variables on each stream in the network. The following conservation of water 
fl-Ux plays an important role 

AeOVeO + AelVel = ^e^e, 6 G /(F). (10) 

the functions A, v and D are defined in the interior of each edge as (for example) 
D{x) = De for a; G e; at internal nodes they are defined as (say) equal to 
zero, and at boundary nodes they are defined as the value of the corresponding 
variable in the node's edge. 

As in assume that mobile individuals follow random trajectories ac- 

cording to a diffusion process on F with drift v and diffusivity D, reflecting 
conditions on the u pstream boundary nodes of the network and absorbing at 



the outlet node (see iFreidlin and Wentzelll . 119931 ). While in the interior of edge 



e, diffusing particles will follow an usual diffusion process with drift and diffu- 
sion coefficients equal to Ue and De respectively. If P{y,x,t), x,y G F, t > 
is the family of transition probability densities for the aforementioned diffusion 
process, then at any internal node e, and as a function of the backwards vari- 
able, P{-,x,t) must be continuous and the total diffusive flux must be equal to 
zero. 

More precisely, for any function / : F — ^ M, and e G /(F), define the total 
flux at e by 

(e) - AeDef'eile) ' Aeoi?eo/eO (0) ^ A^l /?el /el (0) , (11) 
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and consider the infinitesimal operator 

A[f] = DJ'J - vj^ on edge e, e e T, (12) 

with domain 

Dom(^) = {/ G C(f ) n C^{T) : f{<f}) = 0, /'(e) = for all e e C/(r), 

df (13) 
■' (e) = 0, for all e G /(r)}. 



dAD 

For any a; G F, define P{-,x,t) as the solution to the "backwards" equation 

^(.,X,t)^A[Pi;X,t)]. (14) 

As in the one-dimensional case described in Section let cr > and 

define the dispersion kernel by 

IC{y,x)^ / ae-'''P{y,x,t)dt, x, y G F. (15) 
Jo 

Performing a change of time scales t„etD = fJ-t on the analog of equation (IT|), writ- 
ing f := r//i and retaining the symbol t for time, gives the integro-differential 
equation governing a population density u, 

— {x,t) = {f-l)u{x,t) + J K,{y,x)u{y,t)Ay. (16) 

Equation (|16p is to be understood as an evolution equation in the Banach space 
C(F), and thus, one can characterize the stability of the trivial solution it = 
via the spectrum of the operator 

IC[f]:^ J^IC{y,x)fiy)dy. (17) 

More precisely (|Appell et alll2000[) . 

Theorem 1.3. Letojic be the largest eigenvalue of the operator IC given in |J7p . 
The inequality 

f - l-f < 

is sufficient for the exponential stability of the trivial solution u = to U6]) . 

1.3 A related Sturm-Liouville problem 

The dispersion kernel JC is closely related to a Sturm-Liouville problem on F. 
Consider the functions 



p{x) := exp j-y" 



^^^^dyl, ,(.) = #!, (18) 



D{y) '^J' ^ D{x)' 
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where the integral on the definition of p is taken along the unique path connect- 
ing the root 4> and the point x £ F. 

Define the following Sturm-Liouville operator, 

- ~{Pefe)' + qefe On edge e, e e r, (19) 

with domain 

Bom{C) = 1/ e C(f) nC2(r); = ^r all e e /(r)l . (20) 

L qad ) 

The specification of "hydrological" boundary conditions is made considering 
functions on the following class: 

Bh = BH{r) = {f:r^ M; /(</.) = 0, /'(e) = for aU e e UiV)} (21) 

The interest here is mostly on functions that belong to Dom(£)ni3/f = Dom(^). 
The next theorem describes how IC and C are related. 

Theorem 1.4. Let C and IC be defined as in m9\) and |J7[ ) respectively. Then 
u! is an eigenvalue of IC, if and only if v = is a q-eigenvalue of L restricted 
to Bh, namely, there exists a solution u to the following problem 

u e Dom{C) n Bh, C[u] = vqu. (22) 

General theory about the Sturm-Liouville problem on graphs, including the 
existence and representation of solutions. Green's functions, and information 
about eigenvalues and eigenf unctions , can be found in recent l iterat ure. In par- 
ticula r, combining general results by IPokornvi and Prvadievl ( 2004 ) and [Below 
(|l988f) one gets 



Theorem 1.5. All eigenvalues of problem 112^) are real, bounded below by one, 
and form a unbounded discrete set. 

And its immediate corollary. 

Corollary 1.6. 1. If r > fi then the population persists. 

2. Let i'i(r) be the smallest eigenvalue of problem Ii22fl . Then 

rc«t:=M(l tft) (23) 



^i(r) 

is the "critical reproduction rate". For r < rcrit the population will face 
imminent extinction. 
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1.4 Conditions for imminent extinction 



Calculation of the smallest eigenvalue i^i{T) to problem ([^ reduces to finding 
the roots of a determinant, and can be performed for small networks. In gen- 
eral, a more useful approach is to compute bounds for the smallest eigenvalue 
that will lead to easily testable sufficient conditions for imminent extinction or 
persistence. 

Variational techniques for eigenvalue bracketing have b een long used for 
Sturm-Liouville problems (see for example Weinbergeil ( 1974I )). and in particular 



for problems on networks by ICurrie and WatsonI ( 2005 ) . Here we extend and 
adapt the latter work to our particular form of operator C 

Our first result relates the eigenvalues of problem (|22|) to those of a similar 
problem on sub-networks of F. First some notation. For any edge e £ F, denote 
by Fe the tree that has e as its root edge. Namely F = F,,, and if e is a "leaf 
edge" then Fg — e. 

Theorem 1.7. Let e €T, and define functions p'^'^\q^''^ and the operator C^'^^ 
by il8\) . U9\) and I120\) for the subtree Fg. Let z^i(Fe) be the smallest eigenvalue 
of the problem C^"^^ [u\ = lyq^'^^u, u e Dom{C'^''^)nBH{Te). Then z^i(F) =^ z/i(Fe). 

In particular, for an upstream sub-network Fg of F, the respective critical 
reproductive rates satisfy rcrit(F) ^ ?'crit(Fe). This implies that any restriction 
of the population to an upstream sub-habitat can only icrease the minimum 
reproductive rate needed for persistence. Also, the following estimate follows 
from theorem (jl.ip . 

uiiT) < minjl + ^QePe + e is a leaf edge | , (24) 

where Pg := and Qe := are the non-dimensional quantities ([S]) defined 
now on each edge. 

Lower bounds for i^i(F) can be found relating problem (1^^ to the similar 
spectral problem with Dirichlet (absorbing) boundary conditions on all bound- 
ary nodes of F. Define the class of functions 

Bd = S£,(F) = {/ : F R; /(e) = for ah e G dV} (25) 

and consider the spectral problem 

M G Dom(/:) n Be, C[u]^f]qu. (26) 

Much more is known about the properties of the solutions to the spectral 
problem (|26l) than those of problem ((22|) . In particular the set of eigenval- 
ues of is real, discrete, unbounded, and the smallest eigenvalue 771(F) 
has an associated eigenfunction that is st rictly positive in the interior of F 
( Pokornvi and Prvadievl 120041: iBelowl Il988t ). 



The spectral problems (H^J and (|26p with hydrological and Dirichlet bound- 
ary conditions can be related using variational techniques, and in particular one 
gets 

z.i(F) <,7i(F), 
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which doesn't prove to be very useful, since a calculation of rji (F) is as difficult 
as that of fi (F). However, the next key lemma offers a different relation between 
the two problems. 

Lemma 1.8. There exists a network F containing T, and differing from F only 
in the length of leaf edges, such that '7i(F) ^ i^i(F). 



As in lPokornvi and Prvadiev lower bounds for ?7i(F) are possible via 



considerations about the oscillation of eigenfunctions for problem (I26p on F. 
Our main result follows from such techniques. 

Theorem 1.9. The smallest eigenvalue i^i{T) for problem i2^) satisfies 

viiT) > min — ^ + 1, for all e e F. (27) 
And in particular, in view of ^2^, 

< < mm 1 - — — - (28) 



eer 4 + PeQe ^J. e is leaf (Q2 + TT^) + 4Q 

2 Advection-diffusion dynamics on graphs 

For ease of notation, assume that the habitat consists of only three channels: 
two smaller channels converge to form another one. The mathematical model for 
this small river network is F, a graph containg three edges: r, (0) and (1). And 
three nodes: the point where the streams converge r, two upstream boundary 
nodes (0), (1), and the outlet cf). Each channel is in reality a three-dimensional 
domain with coordinates (x, y, z) where the longitudinal variable takes values 
between [0, le], x = corresponding to the downstream end of the channel, and 
le being the length of the channel. 

In full generality, organisms dispersing in a channel e can be quantified by a 
volumetric concentration Ce(a;, y, z, t) given by the number of mobile individuals 
in a representative unit volume of water. Let We{x,y,z,t) denote the water 
velocity field. Assuming that the mechanics of mobile individuals follow a linear 
advection-diffusion process, then there exists a positive-definite diffusion tensor 
Tie{x,y,z,t) such that 

(JC 

^ = V- (DeVCe) - V- (VeCe). (29) 

Several simplifications are in order. First, assume that Ve — (— fe,0,0) 
where Ue > is constant, and the diffusivity tensor is the constant matrix 
De = Delsxs- The transversal variables y, z vary over the cross-sectional area 
of the channel which we will assume constant Af^{x) = A^. The concentration 
of individuals per unit length in edge e, 

u{x,t) := / Ce{x,y,z,t)dydz (30) 
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Figure 3: Schematic representation of a junction of two channels. The corre- 
sponding model is a graph with three edges: r, (0) and (1) connected at node 
r. 



satisfies the usual one-dimensional advection-diffusion equation 

X e e, t>0. 



- — 

dt dx^ 



+ Ve 



dx 



(31) 



The conditions on u at the intersection of the channels are now derived. In 
the interior of each channel, the continuity of the concentration c is implied 
by equation (j29p . The junction of the channels is a three-dimensional volume 
depicted by the shaded area in Figure ([3]). Assuming that c is constant there, 
yields the jump condtion of u at node r. 



Ur{lr,t) _ U(o)(0,i) _ W(i)(0,i) 



A 



(0) 



Ail 



(32) 



Two variables are conserved at the junction: the water fiow, and the flux of 
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mobile individuals. The former condition yields 



A 



(33) 



Also, if storage and other local population changes are neglected, the total 
amount of individuals entering the junction volume must be equal to the amount 
leaving it. Therefore, 



e=(0),(l) 



De^ 1- VeCe 

OX 



LJip ~\- Uf Cf 

ox 



(34) 



which in view of ((33|) and (l32t yields 



du 



(0) 



dx 



(35) 



This may be viewed as an interface condition generalizing notions of "one - 
dimensional skew diff usion" to a tree graph; see lAppuhamillage et all (|201ll ): 
Ramirej ( To appear. ) , and references therein for a discussion of surprising con- 
sequences of these and other such interface conditions on the dispersion of par- 
ticles across an interface in one-dimension. 

Now we address boundary conditions for equation (j3ip on the small river 
network F. First, it is assumed that any mobile organism that reaches the 
outlet (f) leaves the habitat forever. In the diffusion process language, this is 
referred to as an "absorbing boundary" and its formalized by imposing 



U{4),t) = UriO,t) = 0, t>0. 



(36) 



The upstream boundaries of edges (0) and (1) correspond to the headwaters 
of the network. Its assumed that no organism crosses those points, namely the 
total flux is equal to zero. 



A^ 



D. 



dx 



0, e=(0),(l). 



(37) 



Consider now a general network F. Denote by (^*, Dom(^*)) the "forward" 
operator defined by (pij) acting on those functions in C^(F) that satisfy (15^ for 
all nodes in /(F), at 4>, and (|57)) for nodes in J7(F). The formal adjoint 
of (^*,Dom(^*)) is the operator {A,T)om{A)) defined in ([H]). Moreover, in 
Freidlin and Sheul ()2000l ) it is shown that A is the infinitesimal generator of a 
strongly continuous semigroup of linear operators on C{T) corresponding to a 
conservative Markov process X = {X{t) : t ^ 0} with continuous sample paths. 
Namely, if P{y, x,t), x,y € T, t > 0, denotes the family of transition probability 
densities defined as the solution of (fT4| . then the semigroup 



uiy,t)^Tt[h]iy) 



h{x)P{y,x,t)dx, i > 0, 



(38) 
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is the solution to 

du 

— =A[u], u{y,0) = h{y). 



As in the classical case (see iBhattacharva and Wavniire , Il990l) . if one defines 
the adjoint semigroup by 

T:[f]^ Jj{y)P{y,x,t)dx, (39) 

then ^*[T(*[/]] = /. Namely, for a suitable initial condition / ^ with 
/p/dx = 1, Tj*[/] given in ([5^ gives the evolution of the population densi- 
tiy. Moreover, the family of transition probability densities has the following 
properties: P{-,x,t) G Doui{A), P{y,-,t) G Doui{A*) for all x,y €T, t > 0. As 
Py := p{y, ., .) solves ^ = A* [P], limt^o P{y, t) = 5y for ah y G T, then for 
^ C r, P{y, X, t) Ax is equal the fraction of individuals that where initially 
at y and occupy A at time t, namely is equal to P(X(t) G A|X(0) = y). 

Let T be an exponential random variable with mean l/u, and independent 
of the process X. Then 



/C(y, x) := P(X(r) G <lx\Xo = y) - / ae-"'P(y, x, t) di (40) 

Jo 

is the transition probability density of a jump process on F starting at y. 

The definition (j40|) for the dispersion kernel /C can be viewed as a time 
homogenization of the "microscopical" dynamics described by the transition 
probabilities P. Namely, the lC[u] term in equation describes the dispersion 
of organisms at time scales comparable to those of population growth. At this 
scale, the continuous motion individuals in the water column scales-up to a 
Markov jump process on F with infinitesimal generator given by the operator 
on the right hand side of (fT6|) . 

3 Analysis of the dispersion kernel 

Stability results for equation ()16p will follow from relating the dispersion kernel 
IC to the Sturm-Liouville operator C defined in (fT9|) with p and q as in (fT8|) . 
By the definition (IT5|) of /C(y, x) and the strong continuity of of the semigroup 
{Tt-.t-^ 0} it follows that for / G C(f ), 

J K:{-,x) f{x)dx eDomiA), {cr - A) j 1C{- , x) f {x) Ax = cr/. (41) 

Defining the operator 

C:=^{a-A) (42) 
one gets the familiar form of the Sturm-Luiville operator (|19l) on each edge. 
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In iPokornvi and Borovskikhl (|2004l ) it is shown that for very general bound- 
ary and internal node conditions, one can find solutions to Sturm-Liouville prob- 
lems in geometric graphs. In particular, if / G C(f ) then, 



IC{y,x)q{x)f{x) dx 



- fiy)- 



Referring back to the "hydrological conditions" Bh in (HU, one therefore has 
that the solution u G Dom(£) n Bh to C[u] = f can be found via the Green's 
function 

Giy,x):=qix)IC{y,x) (43) 

by making u{y) = J^G{y,x)f(x)dx. More specifically, we have the following 
theorem from which theorem (|1.4p follows as a corollary. 

Theorem 3.1. Let IC{x,y) be as in if75)) . IC[f]{x) — IC{y , x) f {y) dy , and £ 
defined on Dom{C) as in il9\) . i20\) . 



1. /C(-, y) e Dom{A*), IC{y, •) G Dom{C) DBh = Dom{A). 

2. The kernel JC satisfies the following symmetry condition 

p{x)A{x)K.{x,y) =p{y)A{y)lC{y,x). (44) 

3. Let f G C(r), then the function u :~ IC[f] is such that -^u G Dom{C) D 
Bh, and C[j^u] = -^f. 

4- u & C(f) satisfies IC[u] ~ uiu if and only if v := -^u satisfies C[v] — ^v. 

Proof. The first assertion follows from the definition of IC{y,x). For {ii) note 
first that C is self-adjoint with respect to dAD in Dom(£) n Bh- Namely, 

J vC[u]dAD^ J uC[v]dAD, w, w G Dom(£) n B//. (45) 

Let f,g & C^(r), and consider the solutions u and v to C[u\ = f and C[v\ — g 
respectively. Using (HSI) yields 



G{x,y)g{y)f{x)A{x)D{x)dxdy = / / G{y,x)g{y)f{x)A{y)D{y)dxdy. 
T Jv JrJr 

Since f,g are arbitrary, the definition gives the desired result. For (Hi), 
use PSI) in M = IC[f] to get, 

"(a;) , f{y) 



I ^(y'2;)-7^^Y^7-^dx = 
Jr p(,x)A{x) 



p{x)A{x) p{y)A{y) 

and use (I43[) to arrive at an expression involving G. The statement in {iv) 
follows directly from {Hi). □ 
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4 Variational approach 



Consider the measure dAD on F with density dAD = A^De dx on edge e, 
and let L\ ^ (e) denote the space of square-integrable functions in (0, le) with 
respect to the measure A^Dedx. Also, let H^{e) denote the Sobolev space 
of functions such that themselves and their first generalized derivative lie in 
L\jj{e). The analog spaces on all of F are given as direct sums: 

i^D(r):=0il^,(e), H\T):=^H\e) 
eer eer 

The real inner product in L\jj is denoted by (•, ■)ad- 

The operator C in can be extended to L\jj(T) by replacing derivatives 
by generalized derivatives where necessary, an its domain is 

Dom(£) = |/ e C(f ) n H\r); pu' e H\r), -^(e) = 0, e e /(r)| . (46) 

We are interested in solving the following spectral problem 

u e Dom(/:) (IBh, jO-[u] = uqu, (47) 

and in particular, finding the smallest number i^i(F) for which a nontrivial 
solution u of ((i7|) exists. 

Extending the results of Currie and Watson ()20Q5l ). one can reformulate 



problem (|22|) in variational form. Define the associated bilinear form 

J-{u,v) ~ J pu v' + quv dAD , w, u e Dom(J^). (48) 

The domain of J- is given by 

Dom(J") = {ue H\T) n C(f ); u{cj}) = 0}. (49) 

Theorem 4.1. 1. A function u belongs to Dom{C) H Bh and solves C[u\ = 
vqu if and only if u Q Dom{J-) and J-{u,v) — i'{qu,v)AD for all v G 
Dom[F). 

2. The smallest q-eigenvalue of C is 

.,(F). inf (50) 
veDom(j^) [qv, v)ad 

Proof. Let u G Dom(J^) with J-(u,v) = iy{qu,v)AD for all v G Dom(7^). Con- 
sider C(j"(F) — ©egrC^(e), the space of smooth functions on each edge with zero 
values at each node of F. Clearly, C^(F) C Dom(7^) and is dense in L^^(F). 
Consider J^{u,(p) for (p G C^(F). The extension of this functional to an oper- 
ator in L\jj(T) is J^(u, •) — —{pu')' + qu = vqu. In particular, pu' G ff^(F), 
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and since p is smooth, u' can be extended continuously to nodes in U(T). For 
arbitrary v G Dom(J^), integration by parts yields 



T{u, v) = {-{pu'y + qu, v)aD + y^^PeVeU'^AeDc 



eer 



Since the the term on the left hand side and the first term on the right both 
equal v{qu,v)AD, the summation on the right must be equal to zero for all v G 
Dom(J"). This yields -^(e) = for all e e /(F), and u'(e) = for e G U{r). 
The converse statement in (i) follows by integration by parts. Part (ii) follows 
from standard arguments n oting that (|50p imp lies the positive-definitness of the 
form •) — {q-, ■)ad (see Weinbergen (1974), page 38). 

□ 

Recall the definition of a subtree Fg C F with root edge e given just before 
theorem (|1.7p . Let cj)^ be the root node of Fg and define p^'^\q^^'> as in (fTSi) 
where integration is taken with </>g instead </> as lower limit. The operator L'^'^^ 
is then defined via (|19p and ([20| on the graph Fg. 



Proof, of theorem flTTp . Let ux{Te) and u'^"^ G Dom(£('=)) n ^^(Fe) be the 
smallest eigenvalue, and a corresponding eigenfunction of problem (j22l) on Fg. 
Denote by J^'^^ the analogous bilinear form of (|48)) and (|49| on the subtree Fg. 
By theorem (|4.1I) , 

= ^^i(Fe)(q(^)w(^),w(^))^,,(r,), 

where (•, ■)ad(t^) is the inner product on L\^{Te)- Both sides of last equation 
can be multiplied by the constant p{(j}^), and the function u^*^^ extended to all 
F via u :~ Uelr^ G Dom(7^), yielding F{u,u) — viiV e){qu,u)AD, therefore 
r/i(Fe) ^ z^i(F). □ 

Recall the spectral problem with Dirichlet boundary conditions introduced 
in (|26|) . As a last application of the variatio nal formulation, we now g ive an 
eigenvalue bracketing result similar to that of ICurrie and Watson ( 20051 ) . 

Proposition 4.2. Lef 771(F) he the smallest eigenvalue of problem !i26\}. Then 
1.1(F) < 771(F). 

Proof. Consider the following bilinear form Td{u, v) — T{u, v) for u,v G T)om{J-£)) 
Dom(J^) ri-Bij. An argument similar to that of theorem (14. ip characterizes 771(F) 
in terms of Td- In particular, 

771(F) mf < mf = 1.1(F) 

veDomiTo} [qV,V)AD veBomiJ^) (qV,V)AD 

since Dom(J^£)) C Dom(J^). Clearly, equality cannot hold. □ 
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5 Eigenfunctions and oscillation 

We now compute the particular form of eigenfunctfons to problem (j22p . relate 
them to eigenfunctfons of problem (j26p , and use results in oscillation theory to 
show the lower bound estimate in theorem (|1.9p . 

Let m denote the number of edges of F. There are a total of m + 1 nodes, 
distributed as #/(r) = 221^ internal nodes, #C/(r) — ^^^^^ upstream bound- 
ary nodes, and one root node. The continuity, flux-matching and boundary 
conditions in Dom(>C) n Bh correspond to a total of 

3#/(r) + #C/(r) + 1 = 2m 

linear functionals, say {ipi] i — l,...,2m}. Namely, a function / belongs to 
Dom(£) n Bh if and only if ipi[f] — for all i. Also, if e is the i-th edge, every 
solution to {—pu'^y + qUe = lyque on e and extended as zero to all of F, is a 
linear combination of 

/i2^-i(a;;t/) :=e"="lle(x), h2^{x■,ly) ef^'^'leix), z = l,2,...,m (51) 

where 

^ Ve + ^vj - - 1) ^ Ve - v/^;g - ADeaji^ - 1) 

■ 2De ' Pe • 2De ' K ) 

Consider the matrix A with entries 

= i,j = l,...,2m. (53) 

The function det(A(i')) is analytic in ly. Moreover, v is an eigenvalue for problem 
(1221) ifandonlyif z^isarootofdet(A(j/)) = 0, and C = (Ci,...,C2m) G Ker(A) 
is such that 

2m 

u{x; v) = Cjhj{x; v) (54) 

is nonzero. In this case, u(-\ v ) is an eigenfunction corresponding to v (see 
Pokornvi and Borovskikh (I2004D ). 



Lemma (|1.8p follows from analyzing the particular form of u{x] v) at leaf 
edges. Fro the proof, we need some additional notation. First, denote by 



Ve , , . yJ\vl-ADe<j{y-l)\ . 

ae:=— , he{u):= — , e £ F. (55) 

Also, recall the proof of theorem ll.il and define 

v*{e):=l + ^ = l + \QePe. (56) 
Then he{i'*{e)) — and cte = + be, /3e = ae — be for v ^ i^*(e); ag = Ce + i6e, 

Pe — O-e — ibe for V > V*{&)- 
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Proof, of Lemma il.S]) . Let e be a leaf edge. A solution u to problem (|22p 
is given on e by Ue{x\v) — Be"''^ + Ce'^"^ for some constants B,C. Using 
u'^{le) — yields C = —B^e^''^"'~^'\ In particular Ue{x;v) — for some 

a; > 0, if and only if ^ = F^{be{le — x)) where = tanh for ^ ^(e-)* and 
— tan for v > v*{e). Since Cg, 6e > for all v, and > for G [1, z^*(e)], 
then there exists le > le such that Ue{le] v) = 0. 

Let F be the network obtained by extending all leaf edges of F to their 
corresponding l^- Define q, and the operator £ on F via (ITS)) and (IT^ . Let 
u be an eigenfunction corresponding to viiT) and u its extension to F. Then 
u e Boif), and = vi{T)qu. Proposition dOJ gives 771(f) < vi{T). □ 

We now need to introduce some language of oscillation theory for eigen- 
functions on networks. To fix ideas, recall that under very general conditions, 
the eigenfunction of the /c-th eigenvalue for the one-dimensional Sturm-Liouville 
problem Cu — vqu on [0, 1] with u(0) = u{l) = 0, has exactly k — 1 zeros in 
(0,Z). For the analog y of this theory to networks we follow the notation and 
essential results from IPokornvi and Borovskikhl ( 2004 ) adapted to the case of 
tree networks. 

Definition 5.1. An S-zone of a function f : T K in C(F) is a subset 5 C F 
such that \f{x)\ > for all x in the interior of S, and f — on dS. The pair 
(£, v) is said to be non- oscillatory on F if no solution u G Dom{C) to C[u\ — vqu 
has an S-zone in F. In particular, if (£, v) is non-oscillatory on F there cannot 
be a sign-constant solution to problem Ii26\) . 

Theorem 5.2. Let iy*(F) = miueer i^*(e). If v < iy*{T), then (£, iy*(F)) is 
non- oscillatory. 

Proof. By theorem 4.2 in Pokornvi and Borovskikhl (2004), it suffices to show 



that for all v < there exists an strictly positive solution u G Dom(£) to 

£[u] = vqu. Let 1 < v < i^*(F), then < tte < /3e are real numbers for all e G F. 
In particular if B^ > \Ce\ then Ue{x;v) = iJee""^ + Cge^^^ > for all x > 0. 
Assume e is not a leaf edge, and has children edges eO and el. It suffices to 
show that if Be > |Ce|, there exist choices of Beo, Ceo, Bei,Cei such that u{-, v) 
remains positive and satisfies the continuity and matching conditions at node 
e. Using the continuity of u at e and the conservation of water equation (jlOp . 
gives (e) = if and only if 

6e0(Se0 - Geo) + bel{Bel - Cei) = 60(^^6"='= - Cecl^^'^). 

Sincere > |Ce|, then if Bee""'" -Cee'^='= > 0. Let c^o G (0,1), Cei :=l-Ceo 
and choose 

Bei = ^ [u{e] v) + "^ei^^ , Cei = "(e; v) - Bei 

Then Ue[le\ = Ueo{Q; v) = Uei(0; v) = u{e; v) and (e) = 0. Moreover 

1 



\Ce., 2 



u{e; v) - Cei — 

Oei 



< Be 
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and hence Uei{-, v) is positive, i = 0, 1. 



□ 



The proof of theorem (|1.9p follows from the following characterization of 
the smallest eigenvalue for problem (126|): 771 (F) is th e suprcmum of real 77 for 
which {C,rf) is non-oscillatory (Pokornvi and Borovsk ikh. 2004, theorem 5.2.2). 
Since v*{r) = ;/*(r), then (£, i/*(r)) is non-oscillatory on F, and therefore 
v*{T) < ryi(f). Lemma gives v*{T) < uiiT). 



6 Discussion 



Corollary (|1.2p and theorem (|1.9|) give necessary conditions for the imminent 
extinction of populations in river networks under advection-diffusion. Only a 
very particular dispersion mechanism is considered: mobile individuals follow 
realizations of an advection-diffusion process inside F for an exponential time of 
mean i. The boundary conditions are also fixed to be absorbent at the outlet, 
and reflecting at upstream boundary points. 

Although the stability properties of equation ([T5)) may be easily general- 
ized to dispersion kernels with different sets of boundary conditions (see prob- 
lem (j26l) for instance), the methods used here depend heavily on the relation 
between /C and the Sturm-Liouville operator £, and in particular to the flux 
condition = at internal nodes. Generalization of the persistence condi- 
tions to other family of dispersion kernels does not seem feasible. On the other 
hand, the particular form of JC used here has the conceptual advantage of being 
derived from microscopical mechanistic considerations, namely those based on 
the th e ory of diffusion processes on graphs developed bv iFreidlin and Wentzell 
(|l993h : lFreidlin and SheJ (|200flf ). 



The ana lysis of the one-dimen sional case of section ()l.ip confirms the results 
reported on iLutscher et a 1 (|2005h even when boundary conditions are imposed. 



Namely, upstream dispersal events in linear advection-diffusion paths, are likely 
enough to ensure the existence of persistence scenarios, even for high values of 
the water velocity. 

Theorem (II. 7p hints at the benefits that living on a network brings to or- 
ganisms susceptible of being washed out of it. In particular, it shows that the 
population can persist in a river network if upstream refuges are in place, namely 
subgraphs Fg of F where rcrit(re) is low. The detailed effects of the topology of 
F on i^i (F) are much harder to track. One reason for this is that the physical 
variables involved in the dispersal kernel are also dependent on this topology, 
e.g. the conservation of water equation (|10|) . 



The bounds for in (|24p might be useful because of their simplicity, but 
are rather crude. Further work is needed to unravel the relationship between the 
smallest eigenvalue (F) and the physical variables on subgraphs and individual 
edges of F. Also, river networks in nature exhibit multiple and well-documented 
scaling properties between channel indexing schemes (e.g. Horton order) and 
phy sical variables like channel length, mean velocity, and cross-sectional area 
(see lRodriguez-Iturbe and Rinaldo. . .l996il . Deriving the consequences that such 
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properties have on the estimation of extinction/persistence conditions will prove 
much useful. 
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